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Abstract. We have analyzed a sample of galaxies be- 
longing to three clusters: Coma, Abell 85, and Abell 496 
(real galaxies) and a sample of simulated elliptical galaxies 
formed in a hierarchical merging scheme [virtual galaxies). 
We use the Sersic law to describe their light profile. The 
specific entropy (Boltzmann-Gibbs definition) is then cal- 
culated supposing that the galaxies behave as spherical, 
isotropic, one-component systems. We find that, to a good 
approximation (~ 10%), both real and virtual galaxies 
have an almost unique specific entropy. Within this ap- 
proximation the galaxies are distributed in a thin plane 
in the space defined by the three Sersic law parameters, 
which we call the Entropic Plane. A further analysis shows 
that both real and virtual galaxies are in fact located on 
a thin line, therefore indicating the existence of another - 
and yet unknown - physical property, besides the unique- 
ness of the specific entropy. 

A more careful examination of the virtual galaxies 
sample indicates a very small increase of their specific en- 
tropy with merging generation. In a hierarchical scenario, 
this implies a correlation between the specific entropy and 
the total mass, which is indeed seen in our data. The scat- 
ter and tilt of the Entropic Line, defined by Lima Neto 
et al. (1999a), are reduced when this correlation is taken 
into account. Although one cannot distinguish between 
various generations for real galaxies, the distribution of 
their specific entropy is similar to that in the virtual sam- 
ple, suggesting that hierarchical merging processes could 
be an important mechanism in the building of elliptical 
galaxies. 

Key words: Gravitation; Methods: N-body simulations; 
Galaxies: clusters: individual: Coma, Abell 85, Abell 496; 
Galaxies: elliptical and lenticular, cD; Galaxies: evolution; 
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1. Introduction 

Elliptical galaxies present a striking regularity in their 
global luminosity distributions. Within a wide range of 
sizes, the light profile of elliptical galaxies can be described 
by a non-homologous generalization of the de Vaucouleurs 
profile, the Sersic law (e.g. Caon et al. 1993; Graham 
& Colless 1997; Prugniel & Simien 1997). 

This regularity may be understood in terms of a relax- 
ation process: elliptical galaxies seem to be in a quasi- 
equilibrium state, implying that they should obey the 
virial theorem. From the second law of thermodynamics, 
a dynamical system in equilibrium is in a maximum en- 
tropy state. Due to their peculiar properties (long range 
unshielded interactions, equivalence of inertial and gravi- 
tational mass, etc.) the thermodynamics of gravitational 
systems present some difficulties, as well explained in aca- 
demic books (e.g. Saslaw 1985). 

For these systems an equilibrium state is never really 
reached. Various dynamical time scales can be defined: the 
natural dynamical time id ~ 1/ \/4:TtGp (where p is the 
mean density of the system) , the violent relaxation time 
scale, <vR, with <vr ~ ^d, related to the phase mixing pro- 
cess which leads to a quasi-equilibrium state, and a large 
secular time scale igoc, which is related to the slow effects 
of two-body gravitational interactions. It is essentially on 
this scale that one can assert that the equilibrium of a 
self-gravitating system is never established. 

For elliptical galaxies, we have ^vR ~ £ ^scc with 
e « {N/\ogNy^ sa 10~^ (where N is the number of 
particles) . Therefore even if the entropy 5 of a galaxy is 
ever growing on the secular time scale, after violent relax- 
ation we have dS/dtyB. ~ £• Stating that the system is in 
a quasi-equilibrium stage is equivalent to saying that the 
entropy is quasi-constant. 
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However, maximizing the entropy results in an isother- 
mal sphere (Lynden-Bell 1967) which is not valuable either 
from the point of view of physics (divergent total mass) 
or from observations (observed density profiles are steeper 
than the isothermal profile; see also White & Narayan 
1987). 

It is important to note that although there are no ex- 
act stationary entropy states for self-gravitating systems 
(that is, no absolute maximum entropy states), lowest 
energy states may exist, as suggested by Wiechen et al. 
(1988). In order to reach such equilibrium states, the sys- 
tem must necessarily undergo a violent relaxation phase, 
be it through a collapse or a merger. However the final 
configuration reached depends, in principle, on how strong 
the violent relaxation phase was. This raises the interest- 
ing question of how these equilibria, based on minimum 
energy, would relate to the final entropy of the system. 

Numerous works have been devoted to the entropy 
problem (see for instance Merritt 1999). In a previous pa- 
per (Lima Neto et al. 1999a, hereafter LGM), a different 
approach has been adopted: instead of trying to obtain the 
final expected configuration by maximizing the entropy, 
LGM admit the existence of a state of quasi-constant en- 
tropy and calculate this entropy by deriving it from the 
observed light (mass) distribution. In order to compare 
objects of different masses, LGM introduced the specific 
entropy s = S/M, that is the entropy normalized by the 
mass. The specific entropy was then calculated by assum- 
ing that the stars obey the equations of state of an ideal 
gas and using the standard thermodynamical definition of 
the entropy. LGM showed that the galaxies of two clus- 
ters and a group had the same value of s and, therefore, 
that one could derive relative distances between these clus- 
ters using the Sersic profile to model the light distribu- 
tion. LGM suggested that the galaxies having an unique 
s could explain distance indicators based on the shape of 
the brightness profile of galaxies, like those proposed by 
Young & Currie (1994, 1995). 

As in LGM, we will describe the light distribution of 
an elliptical galaxy using a Sersic profile: 



with: 



1.3 



E(i?) = EoCxp(-(i?/a)'') 



(1) 



characterized by three primary parameters: z/, the shape 
parameter (independent of cluster distance), a, the scale 
parameter (distance dependent, in arcsec), and Sqj the 
intensity parameter (in erg arcsec"^). 

In contrast with LGM, who used the thermodynami- 
cal definition of the entropy, we will adopt here the micro- 
scopic Boltzmann-Gibbs definition, therefore eliminating 
the assumption based on the equations of state of an ideal 
gas. Assuming that elliptical galaxies are well described by 
the Sersic law, we have derived the specific entropy (see 
details in Appendix I): 



F{v) = +0.21n(i^) - — + 3.9iy 



-1.3 



- 2.7. 



(3) 



s{a, u, So) = 0.5 ln(I]o) + 2.5 ln(a) + F{iy), 



(2) 



Should the specific entropy of galaxies, s(a,i/. So), be a 
constant or, at least, display a small dispersion around 
its mean value, then Equation (^ would define a thin 
surface in the parameter space [So, a, j/], or a plane in the 
space [In(So), ln(a), F(i^)]. The results presented in LGM 
suggest that this is indeed the case. 

In the next section we describe the data used in this 
paper, i.e. the surface brightness of cluster galaxies as well 
as that of simulated galaxies. We also discuss the fitting 
techniques used to derive the Sersic profile parameters ap- 
pearing in Equation (||) ; in Sections 3T and 3^ we look for 
correlations between these parameters; in Section [3.3| we 
show that besides the uniqueness of the specific entropy of 
galaxies, another relation is also observed. This question 
is revisited in the context of the cosmological scenario of 
hierarchical merging galaxy formation in Section ^. We 
argue that the observed variations of the specific entropy 
of galaxies are correlated with their total luminosity (or 
mass). We then show in Section ||, how this correlation 
helps to understand the tilt of the Entropic Line defined 
in LGM and therefore to further refine the profile-shape 
distance indicator of galaxies based on the shape param- 
eter. We discuss our results in the last section. 

2. The data and fitting methods 

Our goal is to fit the surface brightness of elliptical galax- 
ies with two or three parameters, to search for correlations 
between these parameters and to look for underlying phys- 
ical properties. 

2.1. The data 
2.1.1. Real galaxies 

We have used data on galaxies belonging to three clus- 
ters: Coma, Abell 85 and Abell 496. These galaxies were 
selected: 1) visually as having an elliptical shape on our 
CCD images, and 2) spectroscopically as having a red- 
shift within the corresponding cluster range. The photo- 
metric data are described in Lobo et al. (1997), Slezak et 
al. (1998) and Slezak et al. (1999), and the spectroscopic 
data in Biviano et al. (1995), Durret et al. (1998) and 
Durret et al. (1999) for these three clusters respectively. 
We have determined the growth curve of each galaxy us- 
ing the ELLIPSE task of iraf|^. The growth curves were 
determined with and without background subtraction. 

^ IRAF is the Image Analysis and Reduction Facility made 
available to the astronomical community by the National Op- 
tical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy (AURA), 
Inc., under contract with the U.S. National Science Founda- 
tion. 
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The first two of these clusters have already been an- 
alyzed for the same purpose in a previous paper (LGM), 
but with circular apertures. We present here a new anal- 
ysis of these two clusters, together with a third cluster 
(Abell 496), based on elliptical fits adapted to the geom- 
etry of each galaxy. 

2.1.2. Simulated (virtual) galaxies 

We have used the merger remnants described in Capelato 
et al. (1995, 1997). There are three merger generations: (1) 
the end-products of merging King spheres (with varying 
impact parameter, relative energy and angular momen- 
tum), (2) mergers between first-generation mergers, and 
(3) mergers between second-generation and between first 
and second-generation mergers. 

2.2. The fitting method 
2.2.1. Real galaxies 

For each galaxy in our cluster sample, we have obtained 
the growth curve (integrated luminosity within elliptical 
regions of area e = ttA x B, where A and B are the semi- 
major and semi- minor axes). 

The background contribution was determined individ- 
ually for each galaxy by fitting the last four points of the 
growth curve as a function of the surface by a straight line. 
We checked the robustness of this result by also fitting the 
last 5, 3 and 2 points of the growth curve. 

After determining the background contribution, we 
have subtracted it from the growth curve. Then we deter- 
mined the total luminosity, Ltot, using the last points of 
the growth curve, the half-luminosity (or effective) radius, 
i?eff , and the radius containing 99% of the total luminos- 
ity, i?99. 

We have then fit the growth curves (corrected for the 
sky) using the integrated form of the Sersic law: 

, , 2^a2 /2 /rY\ 

L(i?) = --So7(^-,(^-j j (4) 

where the value of R is not a radius but an equivalent 
radius, R = V AB, and 7(c, x) is the standard incomplete 
Gamma function. The luminosity growth curve fits were 
done with a standard least square minimization method, 
using the 'minuit' programme from the CERN software 
library. 

In order to avoid effects due to the seeing, we have used 
only data points from 2.0 arcsec outwards (the seeing was 
FWHM « 0.9 arcsec for Coma and 1.2 arcsec for A85 and 
A496). The fits were done using data points up to Rgg so 
that for all galaxies the same amount of light was used for 
the fits. The results of these fits are given in Tables |[ ^, 
and|. 



2.2.2. Simulated galaxies 

For simulated galaxies we have computed a mean radial 
mass growth curve (that is, integrated mass instead of 
integrated luminosity), as follows. Each galaxy was ran- 
domly projected in a plane and a growth curve was com- 
puted for each projection by simply counting the particles 
inside iso-density ellipses. These projected growth curves 
(500 for each galaxy) were then used to compute a mean 
one. Having determined a growth curve for each simu- 
lated galaxy, we have proceeded in the same way as for 
real galaxies, fitting it to the integrated Sersic profile. 

The effective radius ReS (i.e. the projected radius con- 
taining half of the total light) is given by LGM: 



Rcs — a R*ff ; 




2 4 6 8 10 

RelT(fit) 



Fig. 1. Effective radius directly measured on the simula- 
tion (Roff(dircct) ) versus effective radius calculated from the 
Sersic fit (Rpff^g).) ) . The straight line indicates the equiva- 
lence between both effective radii. 

We notice that the simulated galaxy profiles we ob- 
tained were extremely close to Sersic profiles. This agrees 
with the results obtained by Capelato et al. (1995) us- 
ing a different fitting technique. In Fig. |^ we compare the 
effective radius of simulated galaxies as estimated from re- 
lation (^ with that directly measured on the simulations. 
As it can be seen there is a good agreement between these 
two quantities. 

3. Relations and correlations 

The last columns of Tables ||, ^, and ^ give the values 
of the specific entropy, sq, for the galaxies of our cluster 
sample. Notice that we cannot compare directly the values 
of So since they depend on a, which is distance dependent 
(we use the apparent a given in arcsec). Neither can we 
compare the values of Sq for a real cluster and for the 
simulation, since they have different units for a and Eg. 
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However, we can compare different generations of mergers, 
as well as the relative dispersion around the mean value 
of So for all real and simulated galaxies. 

Table |l| gives the mean values and dispersions of the 
specific entropy for the whole sample of galaxies. As ex- 
pected, the mean specific entropy varies from cluster to 
cluster, reflecting the different distance of each cluster (cf. 
LGM). For each real cluster, the dispersion is about 10% 
around the mean value. 

For virtual galaxies, the specific entropy seems to in- 
crease with the hierarchy, but with a much smaller disper- 
sion, around 5% of their mean values within each gener- 
ation. This small dispersion is reminiscent of the results 
discussed by Capelato et al. (1997), which show that the 
scatter of the Fundamental Plane defined by virtual galax- 
ies is smaller than for the observed ones by a factor of 
about 2. The increasing of sq with the hierarchy of the 
merger will be addressed in Section ^. 
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Table 1. Mean specific entropy statistics. 
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It is the relatively small scatter of the specific entropy 
of galaxies around their mean values which justifies the re- 
sults discussed in LGM, leading to the definition of a mean 
specific entropy plane, defined through Equation (||). We 
will call it the Entropic Plane. However, as we will see 
in the following, there is another relation linking the ob- 
served quantities of galaxies. 

3. 1 . Correlations with real galaxies 

The correlations of the Sersic profile parameters taken two 
by two are displayed in Figures ||, || and ||. Notice that we 
used the "astronomical like" quantity — 2.51ogEo instead 
of Eo. 

These correlations all have the same general aspect as 
those presented in LGM: the three parameters appear well 
correlated two by two. This strongly suggests that galaxies 
are not distributed randomly around their mean Entropic 
Plane but, instead, are distributed along a thin curve in 
this plane. 

3.2. Correlations with virtual galaxies 

A similar analysis applied to virtual galaxies is shown in 
Figure ^ The correlations of the Sersic parameters taken 
two by two are similar to those found for real galaxies. 
However, the scatter for virtual galaxies is much smaller. 



12 I • • 

0.001 0.01 0.1 1 10 

a (arcsec) 

Fig. 2. The correlation between a and Eq for galaxies be- 
longing to the following clusters: (top) Coma; (middle) 
AbeU 85; (bottom) AbeU 496. The fines are those dis- 
cussed in Section and Appendix II. 

3.3. Another relation? 

The results discussed above suggest that we may go one 
step further: since the projections of the galaxies belong- 
ing to the Entropic Plane are simultaneously three thin 
curves, the galaxies must in fact be located on a thin curve 
in the Entropic Plane. 

In order to check this hypothesis, we will look at the 
Entropic Plane edge-on and face-on. This requires a rota- 
tion to a new coordinate system [^,77, C]i defined by: 

e = [ln(a)-51n(Eo)]/V26 

11 = [-51n(a) - In(Eo) + liF{v)]/^/m> 

C = [51n(a)-}-ln(Eo) + 2F(i.)]/V30. (6) 

We apply this rotation to each galaxy and show the 
result in Fig. ^ for virtual galaxies. This figure suggests 
several comments: 

— The value of the specific entropy depends only on C, 
s — \/7.5 C- Indeed, as a first approximation, virtual 
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galaxies all have the same specific entropy, its numeri- 
cal value depending on the choice of units. Deviations 
from this constant value will be addressed in the next 
section. 

— The galaxies are all effectively located on a curve (in 
fact very nearly a straight line): 

L(f,r;) = L(I]o,a,i/) =0. 

Applying the same rotation to the real galaxies data 
gives quite similar results as for the virtual ones, as seen 
in Fig. 0. However: 

— The scatter is much larger; 

— It is not clear whether the relation L{^,ri) — found 
for virtual galaxies can be approximated by a straight 
line as before. 



The curves displayed in Figs. were obtained by 
assuming that the relation L{^,r]) — may be approxi- 
mated by a straight line. Their derivations are given in 
Apppendix II. 



middle: Abell 85; bottom: AbeU 496. 

At present we do not know of any physical explana- 
tion for this relation. In order to have a curve in a three- 
dimensional space, we need two surfaces to intersect. If one 
of them is the Entropic Plane, then the other one must be 
derived from an independent relation, for instance, a scal- 
ing law relating the gravitational potential energy to the 
total mass, or the depth of the potential well to the mean 
potential (as also suggested by Lima Neto et al. 1999b). 
Such an intersection of surfaces is actually a geometri- 
cal interpretation of the scaling relations that govern the 
physics of a given system. This kind of problem is similar, 
for example, to the one encountered when dealing with the 
origin of the TuUy-Fisher relation (e.g. Mo et al. 1998) or 
of the luminosity-temperature relation for X-ray clusters 
of galaxies (e.g. Markevitch 1998). 

4. Is the specific entropy really unique? 

The merging scenario 

As explained above, the virtual elliptical galaxies we are 
using come from the merging of successive generations. In 
such merging processes the energies and masses of the pro- 
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Fig. 5. The correlations between the Sersic profile param- 
eters for the simulated galaxies. The lines have the same 
meaning as for Figs. ^, ^ and ^. 

genitors will be redistributed in order to generate another 
elliptical. In this section we analyze the effect of merging 
on the value of the specific entropy. 

When viewing a movie displaying the simulation of a 
galaxy merger as compared to one displaying a cold self- 
gravitational coUisionless collapse, one is struck by the 
much more violent matter motions (which are actually 
the engine for the violent relaxation process) occurring in 
the collapse simulation. Furthermore, what one observes in 
such movies is, in a sense, quite similar to an observation 
of the real universe, that is, a macroscopic observation. 
This is in contrast with the microscopic description one 
may obtain from the knowledge of detailed evolution of 
the phase space positions of each particle provided by the 
simulation. In other words, the relevant description of the 
system is done with a coarse-grained distribution function. 

In any case, the mixing in phase-space that occurs dur- 
ing violent relaxation is responsible for the increase of the 
coarse-grained entropy of a dynamical system (Tremaine 
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10 f 
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Fig. 6. Virtual galaxies. Right panel: Entropic Plane seen 
face-on. Left panel: Entropic Plane seen edge-on. Notice 
how which is equivalent to s/^/7.5, has a much smaller 
variation vis-a-vis 77 and ^. 



et al. 1986, Merritt 1999). However, it is not clear how 
much the entropy or, perhaps more important, the spe- 
cific entropy, increases during the violent relaxation phase. 
Does the entropy increase depend on the amplitude of the 
time- varying potential induced by the violent matter mo- 
tions occurring during this phase? It could be that the 
increase of entropy is small compared to the total mass of 
the system (note that defining the entropy as in Eq. (|lO|) , 
S has the unit of a mass). In this case, the specific entropy 
increment could be insignificant when compared with the 
change of other quantities like, e.g., the total gravitational 
energy or the total mass of the system (as indeed will be 
shown below). 

4-2. Shift of the specific entropy 

A careful inspection of the specific entropy (which is given 
by VT.S C) both for virtual and real galaxies indeed shows 
that it does vary, although slightly. We have zoomed the 
left panel of Fig. ^ and show the result in Fig. ||: an over- 
all increase of the specific entropy is observed (C is not 
exactly a constant). Three vertical lines have been drawn, 
corresponding to the mean entropy of the three successive 
generations of galaxies. The specific entropy is actually 
different for each generation of galaxies, and seems to in- 
crease by quanta of specific entropy from one generation 
to the next, although the jump of specific entropy is quite 
small (between 10 and 20 %, see also Table |l|). 

One obvious difference between galaxies of different 
generations is the increase of their total mass. In fact, 
even in a given generation there are galaxies with differ- 
ent masses because of slightly different initial orbital pa- 
rameters of the progenitors. Therefore, we have plotted 
the specific entropy as a function of the total light (mass) , 
i.e., the integrated luminosity (mass) given by relation ^ 
extrapolated to infinity: 



Ltot — So r(— ) 



(7) 
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Fig. 7. Right panel: Entropic Plane seen face-on. Left and 
bottom panels: Entropic Plane seen edge-on. Top: Coma, 
middle: Abell 85, bottom: Abell 496. 



The corresponding Fig. ^ shows that the total mass allows 
to discriminate clearly the entropy of the three generations 
of galaxies. Notice that it is not the mass by itself which 
is really responsible for the shift of the specific entropy, 
but rather the merging process. 

Other parameters do not allow such a clear discrimi- 
nation between generations. For instance we have plotted 
the specific entropy as a function of the v parameter for 
virtual galaxies (Fig. |l0|) . For a given value of i', several 
values of the specific entropy are possible, implying that 
the parameter v is not a good discriminant between gen- 
erations. 

In Fig. |ll| we show a plot similar to Fig. ||, but now 
using the real galaxies data; the same overall trend of an 
increase of the specific entropy (C) with rj seems to oc- 
cur with real galaxies. The striking similarities of these 
two figures further reinforces the hypothesis that ellipti- 
cal galaxies have been formed by mergers. 
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Fig. 8. Enlargement of the (-rj view of the Entropic Plane 
for virtual galaxies. Each merger generation is represented 
by a different symbol. The vertical lines show the mean 
values of ( for each generation. 
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Fig. 9. Virtual galaxies. Specific entropy versus total 
mass. Symbols are the same as in Fig. |[ 

Although we do not know to which generation each 
real galaxy belongs, we may push the analogy between 
real and virtual galaxies a little further by searching for 
correlations between the entropy and luminosity. The cor- 
responding plots are shown in Fig. Again, there is a 
correlation between sq and Ltot for real galaxies, although 
with lower signal-to-noise ratio. We now investigate the 
implications of such a correlation in the context of the 
Entropic Line defined by LGM. 

5. The Entropic Line revisited 

LGM have proposed to rewrite Equation as: 

y(So,a) = 0.51n(Eo) + 2.51n(a); 
AM = F{,,); 
Y + X = so. (8) 

Equation (||), relating 1^ as a function of A, is the equa- 
tion of a straight line with a slope — 1 , called the Entropic 
Line in LGM, which is, in fact, the Entropic Plane seen 
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Fig. 10. Virtual galaxies. Specific entropy versus v. Gen- 
eration symbols are the same as in Fig. ||. 

edge-on; in this equation sq is the mean value of the spe- 
cific entropy of galaxies in a given cluster (as in Table |l|). 
Remember that sq is a distance dependent quantity and 
thus that values for different clusters are not directly com- 
parable. 

In their paper, LGM compared the observed data to 
the predictions of Eq. (||), showing that there is a small 
difference between the corresponding slopes (see Table ||, 
col. 3). The fact that the Entropic Line obtained with the 
observed data had a slope different from —1 was referred 
to as the tilt of the Entropic Line. However, as we have 
seen in the last Section, the specific entropy varies with 
total mass (or light) and, as we will see below, this may 
be at the origin of the tilt. 

Table 2. Statistical analysis: Coma, Abell 85 and Abell 
496. Col. 2: mean specific entropy (note that sq is not 
directly comparable for different clusters, see text); col. 
3: slope of the correlation of Y{T,o,a) with X(i/); col. 4: 
slope of the correlation of F(So,a) with Xcor{i^)', col. 5: 
standard deviation of the data residuals relative to the En- 
tropic Line; col. 6: standard deviation of the data residuals 
relative to the corrected Entropic Line. 
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Fig. 11. Rotated Entropic Plane for real galaxies. Upper 
panel: Coma Cluster. Middle panel: Abell 85. Lower panel: 
Abell 496. Notice that the ( axis is enlarged. This figure is 
comparable to Fig. || (keeping in mind the different units). 



We assume the specific entropy so to be a function of 
the total luminosity, Ltot, as given by Eq. ^ 

So = Q;ln(Ltot) + so,o 

the slope a being obtained through an ordinary least 
square mean fitting of the data (shown as lines in Fig. 
p^ ). The intercept so,o may be understood as a corrected 
specific entropy, that is, taking into account the correla- 
tion with Ltot- 



A "corrected" equation can then be written in place of 
equation (||) as: 

Acor(j^) = F{iy) - a \n{Ltot) 
Y + X^or = so,o (9) 

In Figs. ^ and |l^ we display both the original 
LGM Entropic Line and the corrected Entropic Line, to- 
gether with their corresponding residuals, for the clusters 
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Fig. 12. Correlation between the specific entropy and the 
total luminosity (mass) for real galaxies. Top panel: Coma 
Cluster. Middle panel: AbeU 85. Bottom panel: Abell 496. 

of our sample. Notice that the slope of the corrected En- 
tropic Line is close to —1, indicating that the tilt between 
the data and the predictions based on an unique specific 
entropy for galaxies, has diminished (See Table |[ col. 4). 
Moreover, the dispersion around the corrected Entropic 
Line is improved compared to the dispersion around the 
uncorrected one (compare cols. 5 and 6 in Table ||). 

6. Discussion and conclusions 

We have shown in this paper that the observational data 
for galaxies belonging to clusters confirm the hypothesis 
that their specific entropy is, to first order, unique. As a 
consequence, the galaxies tend to stay in a thin plane (the 




10 

X(v) 



15 



20 



25 



Fig. 13. Coma Cluster. Top panel: Entropic Line and cor- 
rected Entropic Line (arbitrarily separated by a constant 
value to facilitate the comparison) . The straight lines have 
a slope of —1, as predicted by Eqs. (^) and (^. Bottom 
panel: residuals. 



Entropic Plane) in the space of the Sersic light profile pa- 
rameters. Moreover, we have also shown that the slight 
observed variations of the specific entropy of galaxies are 
correlated to their total luminosities. Henceforth, by tak- 
ing this correlation into account, we were able to apply a 
correction which resulted in decreasing the scatter of the 
Entropic Plane. 

We have also shown that, besides the Entropic Plane, 
another relation must exist between the Sersic profile pa- 
rameters of galaxies. The intersection of this relationship 
with the Entropic Plane defines a curve in the 3D Sersic 
parameter space which appears very well defined in our 
data. The physical origin of this new relation is unknown, 
although it may be related to the gravitational energy of 
galaxies, as suggested by Lima Neto et al. (1999b). 

The existence of this curve explains the tight correla- 
tions between the Sersic parameters, as noted in LGM. It 
constitutes the theoretical background for the photometri- 
cal distance indicator proposed by various authors (Young 
& Currie 1994, 1995; Binggeh & Jerjen 1998). 

Our study is solely based on quantities extracted from 
photometric data. The dynamics are in fact hidden in the 
shape parameter of the fitting Sersic profiles. Usually, both 
photometric and spectroscopic studies are performed on 
galaxies, leading for instance to the "Fundamental Plane" . 
It is interesting to notice that some authors have shown 
that globular clusters (Bellazzini 1998), galaxies (Guzman 
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Fig. 14. Same as in Fig. P for Abell 85. 




Fig. 15. Same as in Fig. pj for AbcU 496. 



plus spectroscopy) from the Entropic Plane (photometry 
alone) will be addressed in a forthcoming paper. 

Elliptical galaxies are well described by a three pa- 
rameter profile; since two constraints are acting, only one 
parameter is free. Because the total luminosity of galaxies 
correlates with the shape parameter (Prugniel & Simien 
1997, Binggeli & Jerjen 1998), their light profile should be 
completely defined by their total luminosity. Shape and lu- 
minosity are probably correlated as a consequence of the 
relaxation process; this is consistent with the fact that 
high luminosity galaxies have fiat profiles while those of 
dwarfs are peaked. 

We must stress that the calculation of the specific en- 
tropy requires a good choice for the density profile. Indeed, 
besides an intensity and a length scale, which are the ba- 
sic parameters of all the laws used (for instance the de 
Vaucouleurs law), structural parameters are also neces- 
sary in order to account for the shapes, which refiect the 
dynamical properties of the galaxies. 

In this context, we found that the choice of the Sersic 
law is particularly suitable not only for the sake of sim- 
plicity, but also because it is the bulk of the light which is 
affected by the variations of the shape factor, as pointed by 
Graham & CoUess (1997). One can find in the literature 
other profiles depending on a shape factor, for instance 
the /3-model. However, due to the asymptotic behaviour 
of this model, the entropy and mass calculated in this case 
are essentially located at large radial distances. Besides, 
the values of these two quantities are very sensitive to any 
computational cut-off, and, moreover, the observational 
data in the outskirts of clusters are generally too poor to 
guide the calculation. Therefore, the /3-model is not suit- 
able for such calculations. 

We do not claim that the Sersic profile is the ultimate 
profile of relaxed systems, any more than for instance the 
de Vaucouleurs profile. However, while the precise analyt- 
ical expression for a profile may not be fundamental, we 
understand that taking into account the shape parameters 
is indeed important. It would then be interesting to see if 
"Universal" profiles - as for instance the NEW (Navarro 
et al. 1995) or the Hernquist (1990) profiles - also have 
shape parameters, allowing the calculation of the specific 
entropy. 
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et al. 1993) or even galaxy clusters (Fujita & Takahara 
1999) are located on a line in that plane instead of pop- 
ulating the whole area a priori permitted by the natu- 
ral range of variation of the parameters. The question 
of how to derive the "Fundamental Plane" (photometry 
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Appendix I: Entropy and the Sersic profile 

We use the microscopic Boltzmann-Gibbs definition of the 
specific entropy. 



S/M = -^ J fie)lnf{e)drdv 



(10) 



where M is the total mass and we assume that the func- 
tion / (the coarse-grained distribution function) depends 
only on the energy e. Then, /(e) can be determined by 
the Abel inversion of the density profile (cf. Binney & 
Tremaine 1987) and the specific entropy may be computed 
(see formulae below). 

In order to compute s we have adopted the following 
hypotheses: 



— Spherical symmetry; 

— Isotropy of the velocity distribution; 

— M/L{r) = constant. 

From the 2D mass distribution given by the Sersic pro- 
file (using the hypothesis of constant M/L ratio), we have 
derived a semi-analytical approximation for the 3D mass 
distribution obtained by deprojecting the Sersic profile 
(see LGM): 



Pir) = Po 



r\-P 



(D 



exp(-[r/a]^); 



p = 1.0 - 0.60971/ + 0.05563t/^ . 

From the 3D mass distribution we can compute the 
distribution function and thus the specific entropy. The 
computation can be done numerically but it is cumber- 
some. We have therefore found analytical approximations 
for the specific entropy (see also LGM): 

1 5 
s(a, ly, T,a) ^ -lnY,o + -\na + F{v) ; 



F{v) = 0.21n(i/) 



1.3 



-1.3 



2.7. 



If s(a, Eg) = So = constant, then the above equations 
define a surface, the Specific "Entropy Plane" (in the ap- 
propriate variables). 

Finally, we give for completeness an analytical approx- 
imation for the corresponding magnitude for a given set 
of Sersic parameters (a, v, Eg): 



m = — 2.51ogL(i? oo) = — 2.51ogl]o — 5 log a + m* ; 



TO* = -0.304iy- 1.708iy- 



Appendix II: Analytical formulae for the correla- 
tions 

The projection of the relation L{^,r]) — 0, introduced in 
Section [3.3^ leads to one dimensional curves in the planes 
[So, a], [So, J/] and [a,v\. We give below semi-theoretical 
formulae for these relations: 

1. we assume that L(^,ri) ~ may be approximated by 
a straight line: 



r] = A^ + B, 



(11) 



where the constants A and B are obtained through a 
fitting of the data; 

2. ^ and t] depend on a, v and Sg by relations (^); 

3. we have postulated that the specific entropy is unique 
(theoretical aspect) i.e. C = sqI^/YI) (C also given by 
Eq. (I)). 
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Combining (6) and (^ij) allows us to recover analytical 
formulae, which are displayed below. We have superim- 
posed these calculated curves on each of the correspond- 
ing data in Figs. ^ ^ and ^, for the real galaxies and in 
Fig. H for the simulated ones. 

We obtain for the virtual galaxies the formula: 

a ^ exp[1.771 x (8.256 - 2.5 x logEo)] 



-2.5 log So = ^i:^-^^:^- 0.05041nH- 7.5123 



exp[^- 1:111 - 0.0891n(.) + 1.317] 



(12) 

and successively for Coma: 

a = exp[0.951 x (-18.76- 2.5 x logEo)] 



-2.5 log Eo = — -14? - 0.0922 ln(t/)-h 20.42 



exp[^-l:^- 0.0971n(.) + 1.741] 



(13) 



Abell 85: 



a 



exp[0.811 X (-19.48 - 2.5 x log Eq)] 



496 1 487 
-2.5 log Eo = ^rpj- - 0.0763 In(i^) + 20.483 



a = exp[^:^ - 14¥ - 0.094 ln(j.) + 1.238] 

(14) 

and Abell 496: 

a = exp[1.3131 x (-19.80 - 2.5 x logEo)] 



901 2 702 
-2.5 logEo = ^^-±i^- 0.13861n(i.) + 21.9 



exp[^-^-0.1051n(.) + 1.6255] 

(15) 



Appendix III: Growth curve fitting results 
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Table 3. Growth curve fitting results for Coma. Tlie Godwin et al. (1983) number is indicated as well as the name 
from either NGC or IC catalogues, v and a are the shape and scale parameters fit with the Sersic law. — 2.51og(Eo) is 
the intensity parameter and the magnitude is in the V band. i?eff and /ies are calculated from the primary parameters. 
F(y) is calculated with Eq. (||) and the specific entropy, sq, with Eq. (||). 

GMP Name v log a — 2.51ogEo I/Tot iioff Moff so 
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Table |. Continued. 



GMP Name v log a — 2.51ogEo I/Tot -Reft A^off ^k^") 
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Table 4. Same as Table || for Abell 85 data. The identification number is that of Durret et al. (1998). 
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± 


0.13 


-0.02 


± 


0.08 


19.42 


± 


0.18 


17.49 


± 


0.53 


1.66 ± 0.13 


20.58 


± 


0.29 


-0.07 ± 0.46 


-9.14 


236 


0.55 


± 


0.01 


-0.57 


± 


0.02 


16.97 


± 


0.04 


15.67 


± 


0.12 


2.43 ± 0.04 


19.59 


± 


0.06 


3.16 ± 0.09 


-7.94 


238 


0.38 


± 


0.01 


-1.49 


± 


0.01 


16.64 


± 


0.02 


17.17 


± 


0.09 


2.13 ± 0.03 


20.81 


± 


0.04 


6.99 ± 0.07 


-9.24 


242 


0.98 


± 


0.01 


0.89 


± 


0.00 


20.23 


± 


0.00 


13.75 


± 


0.01 


13.41 ± 0.01 


21.38 


± 


0.00 


-0.08 ± 0.01 


-4.28 


243 


0.32 


± 


0.01 


-1.78 


± 


0.03 


15.19 


± 


0.04 


15.17 


± 


0.19 


4.37 


± 


0.12 


20.37 


± 


0.07 


9.68 ± 0.15 


-7.55 


246 


1.75 


± 


0.62 


0.29 


± 


0.05 


22.77 


± 


0.18 


20.02 


± 


0.40 


1.74 


± 


0.76 


23.22 


± 


0.33 


-1.44 ± 0.51 


-10.28 


253 


0.76 


± 


0.06 


-0.25 


± 


0.08 


18.75 


± 


0.17 


17.29 


± 


0.51 


1.68 


± 


0.12 


20.42 


± 


0.25 


1.01 ± 0.40 


-9.07 


263 


0.48 


± 


0.01 


-1.06 


± 


0.01 


17.29 


± 


0.02 


17.67 


± 


0.05 


1.38 


± 


0.01 


20.37 


± 


0.02 


4.23 ± 0.04 


-9.82 


283 


0.64 


± 


0.01 


-0.17 


± 


0.02 


18.90 


± 


0.03 


16.43 


± 


0.11 


3.29 


± 


0.05 


21.01 


± 


0.05 


1.96 ± 0.09 


-7.73 


305 


1.16 


± 


0.06 


0.11 


± 


0.02 


20.09 


± 


0.06 


17.80 


± 


0.15 


1.72 


± 


0.04 


20.98 


± 


0.09 


-0.61 ± 0.13 


-9.23 


316 


1.11 


± 


0.07 


0.28 


± 


0.03 


19.92 


± 


0.08 


16.73 


± 


0.21 


2.70 


± 


0.08 


20.88 


± 


0.12 


-0.48 ± 0.17 


-8.05 


324 


0.34 


± 


0.01 


-1.43 


± 


0.01 


17.92 


± 


0.02 


16.93 


± 


0.07 


5.72 


± 


0.05 


22.71 


± 


0.03 


8.64 ± 0.06 


-7.84 


326 


0.70 


± 


0.04 


-0.36 


± 


0.06 


18.56 


± 


0.11 


17.38 


± 


0.36 


1.62 


± 


0.08 


20.42 


± 


0.17 


1.42 ± 0.28 


-9.19 


329 


0.81 


± 


0.02 


-0.25 


± 


0.02 


19.94 


± 


0.04 


18.69 


± 


0.12 


1.44 


± 


0.02 


21.48 


± 


0.06 


0.72 ± 0.10 


-9.93 


413 


1.07 


± 


0.08 


0.02 


± 


0.04 


19.94 


± 


0.10 


17.96 


± 


0.28 


1.59 


± 


0.06 


20.96 


± 


0.16 


-0.36 ± 0.24 


-9.41 


447 


0.47 


± 


0.01 


-1.00 


± 


0.01 


16.66 


± 


0.01 


16.61 


± 


0.03 


1.76 


± 


0.01 


19.83 


± 


0.01 


4.47 ± 0.03 


-8.98 



16 I. Marquez et al.: Gravo-thermal properties and formation of elliptical galaxies. 

Table 5. Same as Table || for Abell 496 data. The identification number is that of Durret et al. (1999). 



Id. V log a -2.5 log Eo I/Tot -Reff Mcff F{v) so 

(arcsec) (mag/D") (mag) (arcsec) (mag/D") 



207 


1.18 


± 


0.18 


0.03 


± 


0.08 


20.92 


± 


0.19 


19.06 


± 


0.49 


1.40 


± 


0.12 


21.78 


± 


0.28 


-0.65 


± 


0.41 


-10.12 


216 


1.38 


± 


0.08 


0.37 


± 


0.02 


21.27 


± 


0.04 


17.89 


± 


0.11 


2.59 


± 


0.04 


21.95 


± 


0.08 


-1.02 


± 


0.12 


-8.67 


237 


0.35 


± 


0.00 


-1.86 


± 


0.00 


14.36 


± 


0.01 


15.87 


± 


0.02 


1.68 


± 


0.00 


18.99 


± 


0.01 


8.19 


± 


0.01 


-9.13 


243 


1.18 


± 


0.12 


0.15 


± 


0.05 


21.04 


± 


0.11 


18.56 


± 


0.30 


1.86 


± 


0.09 


21.91 


± 


0.18 


-0.65 


± 


0.28 


-9.46 


243 


0.92 


± 


0.12 


-0.02 


± 


0.09 


20.43 


± 


0.21 


18.37 


± 


0.60 


1.86 


± 


0.16 


21.71 


± 


0.32 


0.16 


± 


0.50 


-9.38 


247 


1.02 


± 


0.07 


0.16 


± 


0.04 


21.50 


± 


0.08 


18.75 


± 


0.24 


2.33 


± 


0.08 


22.58 


± 


0.13 


-0.23 


± 


0.22 


-9.20 


254 


0.41 


± 


0.01 


-1.18 


± 


0.07 


15.65 


± 


0.11 


15.25 


± 


0.42 


2.83 


± 


0.16 


19.51 


± 


0.17 


6.18 


± 


0.32 


-7.84 


257 


1.12 


± 


0.05 


0.18 


± 


0.02 


20.31 


± 


0.05 


17.60 


± 


0.14 


2.14 


± 


0.04 


21.25 


± 


0.09 


-0.51 


± 


0.14 


-8.81 


258 


1.10 


± 


0.04 


0.14 


± 


0.02 


20.13 


± 


0.05 


17.58 


± 


0.13 


2.02 


± 


0.03 


21.10 


± 


0.07 


-0.45 


± 


0.11 


-8.89 


259 


0.32 


± 


0.00 


-1.90 


± 


0.03 


14.55 


± 


0.05 


15.00 


± 


0.21 


3.64 


± 


0.11 


19.80 


± 


0.08 


9.89 


± 


0.16 


-7.78 


260 


0.89 


± 


0.06 


-0.09 


± 


0.05 


19.48 


± 


0.12 


17.67 


± 


0.33 


1.68 


± 


0.08 


20.80 


± 


0.17 


0.26 


± 


0.27 


-9.21 


261 


0.94 


± 


0.07 


-0.08 


± 


0.06 


19.03 


± 


0.13 


17.32 


± 


0.35 


1.55 


± 


0.08 


20.27 


± 


0.19 


0.08 


± 


0.28 


-9.16 


262 


0.83 


± 


0.07 


-0.14 


± 


0.07 


19.36 


± 


0.16 


17.59 


± 


0.46 


1.79 


± 


0.12 


20.85 


± 


0.23 


0.60 


± 


0.36 


-9.10 


264 


0.53 


± 


0.01 


-0.67 


± 


0.02 


17.06 


± 


0.04 


16.11 


± 


0.13 


2.16 


± 


0.04 


19.77 


± 


0.06 


3.36 


± 


0.10 


-8.34 


266 


0.83 


± 


0.03 


-0.29 


± 


0.04 


17.97 


± 


0.08 


17.00 


± 


0.23 


1.23 


± 


0.04 


19.44 


± 


0.12 


0.58 


± 


0.17 


-9.39 


267 


0.90 


± 


0.07 


-0.12 


± 


0.06 


18.46 


± 


0.15 


16.84 


± 


0.40 


1.54 


± 


0.09 


19.77 


± 


0.21 


0.24 


± 


0.31 


-8.97 


268 


0.83 


± 


0.03 


-0.23 


± 


0.03 


18.95 


± 


0.07 


17.67 


± 


0.20 


1.41 


± 


0.04 


20.41 


± 


0.10 


0.56 


± 


0.15 


-9.49 


272 


0.91 


± 


0.07 


0.03 


± 


0.06 


21.05 


± 


0.12 


18.71 


± 


0.36 


2.12 


± 


0.11 


22.34 


± 


0.19 


0.20 


± 


0.30 


-9.34 


287 


0.55 


± 


0.01 


-0.54 


± 


0.03 


17.38 


± 


0.05 


15.97 


± 


0.19 


2.52 


± 


0.06 


19.98 


± 


0.08 


3.08 


± 


0.15 


-8.02 


288 


0.68 


± 


0.01 


-0.32 


± 


0.02 


17.75 


± 


0.04 


16.27 


± 


0.14 


1.92 


± 


0.04 


19.69 


± 


0.07 


1.58 


± 


0.11 


-8.44 


291 


0.41 


± 


0.01 


-0.98 


± 


0.03 


16.08 


± 


0.05 


14.81 


± 


0.20 


4.07 


± 


0.11 


19.85 


± 


0.08 


5.97 


± 


0.16 


-7.06 


293 


0.96 


± 


0.03 


0.03 


± 


0.02 


19.95 


± 


0.05 


17.74 


± 


0.14 


1.91 


± 


0.04 


21.13 


± 


0.08 


-0.02 


± 


0.12 


-9.05 


294 


0.47 


± 


0.00 


-0.85 


± 


0.02 


16.07 


± 


0.03 


15.15 


± 


0.10 


2.70 


± 


0.04 


19.30 


± 


0.04 


4.60 


± 


0.08 


-7.69 


295 


0.87 


± 


0.04 


-0.03 


± 


0.04 


19.82 


± 


0.09 


17.64 


± 


0.26 


2.05 


± 


0.07 


21.20 


± 


0.14 


0.39 


± 


0.22 


-8.91 


304 


1.00 


± 


0.04 


0.08 


± 


0.03 


19.75 


± 


0.06 


17.37 


± 


0.17 


2.00 


± 


0.05 


20.87 


± 


0.09 


-0.14 


± 


0.15 


-8.80 


306 


0.78 


± 


0.01 


-0.15 


± 


0.02 


18.43 


± 


0.04 


16.58 


± 


0.12 


1.96 


± 


0.03 


20.03 


± 


0.06 


0.85 


± 


0.09 


-8.51 


309 


0.44 


± 


0.01 


-0.93 


± 


0.08 


16.48 


± 


0.12 


15.41 


± 


0.49 


3.33 


± 


0.22 


20.01 


± 


0.19 


5.37 


± 


0.39 


-7.57 


311 


0.85 


± 


0.03 


0.00 


± 


0.03 


20.37 


± 


0.06 


18.00 


± 


0.20 


2.30 


± 


0.06 


21.81 


± 


0.10 


0.50 


± 


0.17 


-8.90 


313 


0.97 


± 


0.04 


0.10 


± 


0.03 


20.55 


± 


0.06 


17.99 


± 


0.18 


2.22 


± 


0.05 


21.72 


± 


0.10 


-0.05 


± 


0.16 


-8.93 


319 


1.00 


± 


0.13 


-0.03 


± 


0.09 


21.01 


± 


0.20 


19.15 


± 


0.54 


1.57 


± 


0.13 


22.13 


± 


0.29 


-0.16 


± 


0.44 


-9.98 


322 


0.50 


± 


0.03 


-0.70 


± 


0.12 


18.16 


± 


0.19 


16.91 


± 


0.74 


2.80 


± 


0.28 


21.14 


± 


0.31 


4.00 


± 


0.58 


-8.38 


326 


0.46 


± 


0.01 


-0.93 


± 


0.07 


17.17 


± 


0.11 


16.54 


± 


0.42 


2.44 


± 


0.14 


20.48 


± 


0.17 


4.79 


± 


0.33 


-8.50 


331 


0.82 


± 


0.02 


-0.04 


± 


0.02 


18.76 


± 


0.04 


16.49 


± 


0.14 


2.26 


± 


0.04 


20.25 


± 


0.07 


0.61 


± 


0.12 


-8.24 


333 


0.41 


± 


0.00 


-0.86 


± 


0.01 


16.58 


± 


0.02 


14.63 


± 


0.07 


5.66 


± 


0.06 


20.40 


± 


0.03 


6.09 


± 


0.06 


-6.50 



